Two-dimensional randomly-frustrated spin 1/2 Heisenberg model 
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We investigate the properties of 5 = 1/2 Heisenberg clusters with random frustration using exact 
diagonalizations. This is a model for a quantum spin glass. We show that the average ground state 
spin is S oc yN, where -/V is the number of sites. We also calculate the magnetic susceptibility and 
the spin stiffness and low-energy excitations, and discuss these in terms of a semiclassical picture. 



PACS: 75.10.Jm, 75.10.Nr, 75.40.Cx 

Frustrated quantum spin systems continue to attract 
much interest. Such systems show rich and complex be- 
havior, including complex ground states, with or without 
long-range magnetic order, quantum phase transitions at 
T=0, and novel excitations/quasiparticles. Most work 
has focussed on model systems with translation symme- 
try, in which the frustration is geometric in origin, result- 
ing from competing antiferromagnetic exchange interac- 
tions. Typical examples are the triangular lattice [jl]]|] 
and the J\ — J2 model on the square lattice Wc 
refer to such systems as "regularly frustrated". A recent 
review of phase transitions in frustrated antiferromagnets 
H discusses many of these issues, particularly for 3-d sys- 
tems. A general field theoretic approach is discussed in 
recent book 0. 

In contrast to these cases are systems where defects 
and/or random interactions or fields lead to frustration. 
Doty and Fisher || investigated the effects of both ran- 
dom fields and random exchange in the spin- 1/2 XXZ 
chain, and found a quantum phase transition controlled 
by exchange anisotropy. Sandvik ||] considered a square- 
lattice antiferromanget with up to 10% ferromagnetic 
bonds, but was unable to reach the interesting spin-glass 
regime. There have also been a number of studies of ran- 
domly diluted systems, but since these are not, in general, 
frustrated we do not discuss them further. 

The most extreme case of frustration by exchange dis- 
order occurs in spin-glass systems [To[ . There is an enor- 
mous body of theoretical work on models of spin-glasses, 
but the bulk of this is on Ising or classical XY and Heisen- 
berg models. The earliest work on quantum spin-glasses, 
an attempt to apply replica theory to the spin-S quantum 
model with long-ranged Gaussian interactions, is due to 
Bray and Moore fLl] . The same model has been con- 
sidered more recently via an SU(N) generalization Jl2| , 
but the relationship of these studies to the more realis- 
tic short-range interaction case remains unclear. To our 
knowledge the only study of a finite S short-range spin- 
glass model is the work of Nonomura and Ozeki Jlj| , who 
carried out exact diagonalizations on small clusters. 

The present Letter is strongly motivated by Ref. [ p^[ , 
and extends that work in several ways. We consider small 
clusters of S=l/2 spins with Hamiltonian 



where Jij are random exchange interactions. Two models 
for randomness are used: the ± J model with 
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with mean and standard deviation 1. We use exact 
diagonalizations for clusters of N=10,16,18,20 spins with 
periodic boundary conditions. These clusters have the 
symmetry of the square lattice, and hence the results can 
be extrapolated, using finite size scaling, to the infinite 
lattice. We have also carried out some calculations with 
free boundaries, in which case our results may be relevant 
to magnetic nanoparticles. 

The first question to be addressed is the average spin 
of the ground state. Using at least 500 different ran- 
domly selected bond configurations for each cluster size, 
we determine the ground state spin quantum number S 
for each realization, and the mean value and standard de- 
viation. Of course, for the homogeneous antiferromagnet 
S=0 rigorously, but for random interactions we find S val- 
ues from to 6 (for the larger clusters). Figure 1 shows 
histograms for N=18 (periodic boundary conditions and 
2000 samples) for both the ± J and Gaussian cases, to- 
gether with the distribution which would be expected in 
the absence of interactions (the number of combinations 
by which one can add 18 spins 1/2 to get the total spin 
S). The histogram for the Gaussian model is slightly 
shifted to larger S compared to that for the ±J model, 
and the pure statistical distribution is shifted further up. 
However, all in all the distributions are very similar. 

From the statistics of ground state spin quantum num- 
bers we have computed the average value of S 2 . One 
can say that this is an average value of the maximum z- 
projection of the spin squared. This quantity has a more 
straightforward semiclassical meaning than < S 2 > — 
S(S+1). In the thermodynamic limit, N — > 00 these 
two quantities coincide. The value of S 2 is plotted in 
Fig. 2 versus N, the number of spins in the cluster. While 
there is some scatter, due both to limited statistics and 
to the different orientations of clusters, the data is con- 
sistent with a linear dependence of S 2 on N shown by the 
dashed lines. In particular 
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± J model : S 2 w 0.120AT, 
Gaussian model : S 2 w 0.202iV. 



This means that the average spin per site varies as 1/ v N, 
and vanishes in the thermodynamic limit. This N depen- 
dence is exactly what would be expected from classical 
fluctuation theory, for a system of spins with random ori- 
entations, though our system is in the extreme quantum 
limit. If we denote the length of the effective classical spin 
by S eff then < S 2 >= N S 2 ff . This yields S e //=0.35 
for the ± J case and 0.45 for the Gaussian case, compared 
to 1/2 for noninteracting classical spins. The quantity is 
reduced due to quantum fluctuations, and we observe 
that quantum fluctuations are greater in the ± J model 
than in the Gaussian case. 

From the same calculations we obtain the ground state 
energy per spin E /N, averaged over at least 500 inde- 
pendent configurations. These values are given in Ta- 
ble 1. There is little variation between different clusters, 
indicating that these results are close to the thermody- 
namic limit. As expected, the Gaussian model yields 
slightly higher energies. 

We have also computed the average spin-glass order pa- 
rameter in the ground state m 2 g = jjrJ2i j < Si ■ Sj > 2 , 
where the line denotes a configurational average. The 
values of this quantity, averaged over 500 random bond 
configurations, are given in Table 1. Plotting the re- 
sults versus 1/N shows a linear variation, which, when 
extrapolated to the bulk limit, yields m sg — 0.01 for 
both ±J and Gaussian cases. This value agrees with 
the estimate in Ref. ]jl3). The effective spin length 
S e ff can be estimated from the spin-glass order param- 
eter, via m 2 sg — S^ff < co 
S e ff — 0.42, consistent with the estimate above. 

As well as the ground state we have studied properties 
of the low-lying excitations of the finite clusters. We have 
determined numerically the spin quantum number of the 
lowest excited state and the energy gap for at least 500 
realizations of bond configurations for N=10, 16, 18,20 as 
before. We observe that predominantly (typically 97% 
of cases) the spin of the first excitation is Si = Sq ± 1 
where So is the spin of the ground state. In Table 1 we 
present /+ and /_: the fraction of configurations which 
have Si = So + 1 and which have S± = So — 1, f° r both ±«7 
and Gaussian models. We also give the average energy 
gap AE for both models, for those cases where Si = 
So ± 1. The average gap for Si = So + 1 is practically the 
same as that for Si = So — 1. The frequency /_ appears 
systematically less than / + . We give an explanation of 
this below. 

The results for the lowest excitations can be under- 
stood within a semi-classical picture, valid for large N. 
We denote by E and |0 > the energy of the ground state, 
which has total spin So oc V^V, and the corresponding 
wave function with maximum S 2 = So- We now consider 
an external magnetic field B in the z-direction. The sys- 



tem will adopt a new ground state, with energy 
E B = Ea-S B-lxNB 2 1 
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>= Sgjy/3, which gives 



where x is the magnetic susceptibility per site. The state 
|0 >, which will precess about the z-axis with Larmor fre- 
quency u) — B, is a "rigid body" excitation in which the 
spins precess coherently. Hence its energy with respect 
to Eb is simply SE = SoB + \lu 2 . Comparing this with 
(||) gives the effective moment of inertia / = N\- Simi- 
lar considerations can be also applied to the 2D quantum 
antiferromagnet and lead to the known results p4| ]. 

The lowest excitations thus correspond to the states of 
a 2D quantum rotor with energies 

L 2 L 2 

E L -E = — = . (6) 

21 2N X 

where L = 0,±1,±2, ... is the angular momentum. The 
total spin of the corresponding excited state is S = Sq+L. 
Thus the lowest excitations will have Si = So ± 1, as 
observed in the numerical results. Note that eq. (||) 
describes rotations around the z-axis. This is why it is 
different from that for the antiferromagnet |l4j which 
describes 3D rotations. In addition to the requirement 
N ^> 1 for this approach, we also require \L\ <C So ~ yN 
to justify the formula (^). In our numerical calculations 
this inequality is only weakly satisfied, and we believe this 
is the reason for the number of samples with Si = So — 1 
being somewhat less than with Si = So + 1. 

In the above consideration the susceptibility \ corre- 
sponds to the magnetic field directed along the total clus- 
ter spin. We argue that the susceptibility determined in 
this way from the cluster data extrapolates, in the ther- 
modynamic limit, to the total susceptibility and not to 
X± (as in the antiferromagnet). This is because, in the 
N — > oo limit, the B 2 term in (||) will be dominant, and 
the energy will be independent of the spin direction. 

From (||) we also see that the energy gap should behave 
like 1/N. In Fig. 3 we plot the quantity 1/2% = N A~E 
versus 1/N. Although finite size corrections appear to be 
large, a rough extrapolation to N = oo gives for the mag- 
netic susceptibility \ ~ 0.40 (±J model), and \ ~ 0-91 
(Gaussian). These values are about an order of magni- 
tude higher than for the 2D quantum S=l/2 antiferro- 
magnet. Another interesting point is that the finite size 
correction is of the opposite sign compared to that for 
the 2D quantum antiferromagnet. 

The final quantity we have computed is the spin- 
stiffness p s , which is a measure of the rigidity of the 
ground state to a small twist 8, 



SE = T^Ps 



(vey d 2 r, 



(7) 



where Eq is the ground state energy. Non-zero spin- 
stiffness, in the thermodynamic limit, indicates a long- 
range magnetic order in the system. Previous calcula- 
tions of p s for finite antiferromagnetic clusters pa] have 
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used a modified Hamiltonian in which the quantization 
axis in the twist direction is rotated by AO for neighbor- 
ing sites, leading to phase factors in the exchange con- 
stants. This method of calculation assumes that the sys- 
tem has an intrinsic spatial periodicity, so it is applicable 
to ferromagnetic or antiferromagnetic states. However 
we found that it cannot be applied to a random system. 
This is why we use a different approach. As the exact 
diagonalization method yields not only the ground state 
energy but also the ground state wavefunction <3> for each 
cluster, we impose the twist directly on the wave func- 
tion, ^0 — Ue 1 ^, where Ug is the unitary transformation 
which rotates each subsequent spin along a given direc- 
tion by an additional angle A6*. The change of the energy 
under the twist, SE =< ^fg\H\^fg > —E Q is calculated 
explicitly, and comparison with (R) gives the spin stiff- 
ness. To avoid ambiguity we have chosen clusters with 
free boundaries, and with the twist imposed along the x 
direction. Table 2 presents values of p s for various clus- 
ters, averaged over 500 bond configurations. There is a 
weak decrease of p s with size of the cluster, but the scat- 
ter due to different shapes does not allow us to make a 
reliable extrapolation. Therefore we take an average over 
clusters value as an estimate for the stiffness. We have 
also performed similar calculations for antiferromagnetic 
clusters and find that the data is very similar: a small 
scattering and a weak decrease with size. However in this 
case we know the thermodynamic value and hence find 
that our cluster calculations overestimate the stiffness by 
the factor 1.4. Assuming that the finite size scaling factor 
for the random system is the same we come to the fol- 
lowing estimates for the stiffness, p s « 0.15 (±J model), 
and p s w 0.14 (Gaussian model). These values are just 
slightly less than the spin stiffness for the 2D quantum 
antifcrromagnet on the square lattice. 

Having the spin stiffness and the magnetic suscepti- 
bility we find the velocity of the Goldstone spin wave 
c = y/ps/x, (see Ref. |l|@), c « 0.6 (±J model) and 
c w 0.4 (Gaussian model). This should be compared 
with c = 1.67 for the 2D quantum antifcrromagnet on 
the square lattice. The low temperature specific heat will 
be due to the Goldstone spin waves and thus will have a 
T 2 temperature dependence. We propose that this will 
be generic for disordered isotropic 2D quantum spin sys- 
tems, and explains, for example, the data for the S = 3/2 
Kagome spin glass |17|| . This should be contrasted with 
Ising spin glasses where the specific heat is due to the 
localized excitations, and is linear in T. 

In summary, we have investigated aspects of highly 
disordered S=l/2 quantum spin systems at zero temper- 
ature by means of exact diagonalizations of small clusters 
of up to 20 spins. Averaging over disorder, we find: l)the 
total spin quantum number scales as S cx V~N, 2)there 
is evidence of spin-glass order in the ground state, in the 
thermodynamic limit, 3)the lowest excitations and en- 
ergy gap for clusters are consistent with the predictions 



of a semi-classical theory, 4)the spin-stiffness does not 
appear to scale to zero in the thermodynamic limit. Two 
models for the random exchange have been used, a ± J 
model and one with Gaussian exchange. There appears 
to be no significant qualitative difference between these 
two cases. Although we have used square lattice clus- 
ters with nearest-neighbor interactions and with periodic 
boundary conditions, to model the behavior of an infinite 
square lattice, we believe our results are more general. 

There are two ways in which our results could be re- 
lated to real systems. The first is to quantum spin- 
glasses, which have received little attention in the litera- 
ture. However, one needs to be cautious in extrapolating 
from small clusters to systems as complex as spin-glasses. 
An alternative application might be to disordered mag- 
netic nanoparticles. We know of no immediate candi- 
dates. However recent work on Mni2 and Fei3 03 spin 
clusters Q (which are not disordered) suggests that this 
is not impossible. 

Finally we remark on a recent study of a system of N 
randomly interacting fermions fl9|j , in which the authors 
argue that ground states with zero and maximum spin 
will dominate. The evidence from our work is that this 
effect does not occur in spin clusters. 

We gratefully acknowledge discussions and correspon- 
dence with C. Hamer, A. Sandvik, and V. Zelevinsky and 
support from the Australian Research Council. 
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FIG. 1. The ground state spin distributions (N=1S, 2000 
samples) for both the ±J and Gaussian models, together with 
the distribution which would be expected in the absence of in- 
teractions. 
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FIG. 2. The average ground state spin squared versus size 
of the cluster. The triangles correspond to the ± J model, and 
the squares correspond to the Gaussian model. 
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AE multiplied by the size of the 
cluster versus inverse size of the cluster. The triangles cor- 
respond to the ±J model, and the squares correspond to the 
Gaussian model. 



Table I. 





N 


10 


16 


18 


20 


±J 


E /N 


-0.493 


-0.498 


-0.500 


-0.500 


model 


m % 


0.051 


0.0355 


0.0325 


0.0303 




U 


0.552 


0.528 


0.504 


0.564 




f- 


0.422 


0.454 


0.464 


0.418 




AE 


0.234 


0.125 


0.100 


0.091 


Gaussian 


E /N 


-0.467 


-0.471 


-0.472 


-0.473 


model 


m% 


0.050 


0.0358 


0.0332 


0.0306 




f+ 


0.512 


0.484 


0.466 


0.510 




f- 


0.474 


0.504 


0.512 


0.462 




AE 


0.179 


0.078 


0.067 


0.057 



The data for different cluster sizes N and for two types of 
random distribution. Eq/N is the average ground state 
energy per site; m 2 sg is the spin glass order parameter; 
/+i is the fraction of configurations which have spin of 
the first excitation So + 1, where Sq is spin of the ground 
state; f-i is the fraction of configurations which have 
spin of the first excitation Sq — 1; AE is the average 
energy gap. 



Table II. 



L X X Ly 


4x3 


6x3 


3x4 


4x4 


5x4 


±J 

model 


0.220 


0.208 


0.214 


0.200 


0.198 


Gaussian 
model 


0.194 


0.190 


0.190 


0.189 


0.190 



The average spin stiffness p s for different rectangular 
clusters. The twist is imposed along the x direction. 
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